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Abstract 

We study numerically the Cauchy problem for equivariant wave maps from 3+1 Minkowski 
spacetime into the 3-sphere. On the basis of numerical evidence combined with stability 
analysis of self-similar solutions we formulate two conjectures. The first conjecture states 
that singularities which are produced in the evolution of sufficiently large initial data are 
approached in a universal manner given by the profile of a stable self-similar solution. The 
second conjecture states that the codimension- one stable manifold of a self-similar solution 
with exactly one instability determines the threshold of singularity formation for a large 
class of initial data. Our results can be considered as a toy-model for some aspects of the 
critical behavior in formation of black holes. 

1 Introduction 

Let M be a spacetime with metric rj and N be a Riemannian manifold with metric g. The wave 
map U : M — > N is defined as a critical point of the action 

S(U)= [ 9AB d a U A d b U B r l ab dVM. (1) 

J M 

The associated Euler-Lagrange equations constitute the system of semilinear wave equations 

d a d a U A + Ti c {U)d a U B d a U c = 0, (2) 

where T's are the Christoffel symbols of the target metric g. 

The recent surge of interest in wave maps (known in the physics literature as cr-models) stems 
from the fact that they provide an attractive toy-model for more complicated relativistic field 
equations. In particular they share some features with the Einstein equations so understanding 
the problems of global existence and formation of singularities for wave maps may shed some 
light on the analogous, but much more difficult, problems in general relativity. Having this 
analogy in mind we have studied the Cauchy problem for Eq.(2) in the case where the domain 
space is 3 + 1 Minkowski spacetime, M — R 3+1 , and the target space is the 3-sphere, N = S 3 . 
Although our primary motivation was an attempt to get insight into some aspects of critical 
behavior in gravitational collapse, we think that our results are interesting in their own right. 
In this paper we restrict attention to equivariant maps. In polar coordinates on R 3+1 and S 3 
the respective metrics are 

?? = -dt 2 + dr 2 + r 2 duj 2 , g = du 2 + sin 2 (u)dtt 2 , (3) 
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where dus 2 and dfl 2 are the standard metrics on S* 2 . Equivariant maps have the form 



U(t,r,u) = (u(t,r),£l(*>)), (4) 

where Q is a homogeneous harmonic polynomial of degree / > 0. In what follows we consider 
the case 1 = 1, where Q = u (such maps are called corotational) . For a corotational map the 
Cauchy problem for Eq.(||) reduces to the semilinear wave equation 

2 sin(2 M ) 

u u = u rr + -u r -J — (5) 

with initial data 



u(0,r) = #r), Ut(0,r)=V(r). (6) 
The conserved energy associated with solutions of this equation 

1 f°° 

E[u] = - / {r 2 u 2 + r 2 u 2 r + 2 sin 2 w) dr (7) 

is manifestly nonnegative and scales as E[u(x/X)] = XE[u(x)] which means that Eq.(||) is su- 
percritical (like Einstein's equations). It is widely believed that for supercritical equations the 
solutions with sufficiently small initial data exist for all times while large data solutions develop 
singularities in finite time 0. In the case of (||) the global existence for small (in the Sobolev 
space H k with sufficiently large k) data was proved by Kovalyov M and Sideris 0. For large 
data there are no rigorous results, however it is known that there exist smooth data which lead 
to blowup in finite time. An example of such data is due to Shatah who showed that (||) admits 
a self-similar solution u(t,r) = foiq&i) which is perfectly smooth for t < T but breaks down at 
t = T. Turok and Spergel JT3j found this solution in closed form /o = 2arctan(j^-) so in the 
following we will refer to /o as the TS solution. 

On the basis of our numerical simulations we conjecture that the example of blowup given 
by Shatah is generic. By this we mean that there is a large open set of initial data which blow 
up in a finite time T and the asymptotic shape of solutions near the blowup point (r = 0) 
approaches /o( ^37) as t — > T~ . In this sense the blowup can be considered as local convergence 
to the TS solution /q. Actually, our failed efforts to produce a non-self-similar singularity lead 
us to suspect that the blowup is universally self-similar. Note that the self-similarity of blowup 
excludes a concentration of energy at the singularity and suggests that the solutions can be 
continued beyond the blowup time in an almost continuous fashion. We must admit that this 
aspect of singularity formation for wave maps is somewhat disappointing from the standpoint of 
modeling the formation of energy trapping singularities (like black holes). 

Whenever the singularities develop from some but not all data, there arises a natural ques- 
tion of determining the threshold of singularity formation. We investigated this issue using a 
basic technique of evolving various one-parameter families of initial data which interpolate be- 
tween global existence (dispersion) and blowup. Along each such family there exists a point 
(critical initial data) which separates the two regimes. We show that the critical initial data 
blow up in a finite time T and the asymptotic shape of solutions near the blowup point (r = 0) 
approaches /i(yzr), a self-similar solution with one unstable mode. The marginally critical 
data approach /j for intermediate times but eventually the unstable mode becomes dominant 
and ejects the solutions towards dispersion or stable blowup (that is, /o). Thus, we conjecture 
that the codimension-one stable manifold of the solution fx plays the role of the threshold of 
singularity formation for a large set of initial data. 

The threshold behavior in our model is similar to the type II critical behavior in gravitational 
collapse (see || for the recent review) where self-similar solutions of Einstein's equations (contin- 
uous or discrete, depending on a model) sit at the threshold of formation of a black hole. There 
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are also many parallels between our results and the work of Brenner et al |4|] on the chemotaxis 
equation. All that suggests that self-similar behavior at the threshold of singularity formation 
is a common feature for evolutionary partial differential equation. 

The rest of the paper in organized as follows. In the next two sections we discuss some 
special solutions of (0) which are potential candidates for attractors. In Section 2 we analyze 
self-similar solutions and their linear stability. Section 3 is devoted to static solutions. In 
Section 4 we describe the results of numerical simulations and document the numerical evidence 
behind the two conjectures formulated above. 



2 Self-similar solutions 



Note that Eq.(|^) is invariant under dilations: if u(t, r) is a solution, so is u a (t, r) = u(at, ar). It 
is thus natural to look for self-similar solutions of the form 



u(t,r) = f 



T-t 



(8) 



where T is a positive constant. Substituting the ansatz (||) into (|J) we obtain the ordinary 
differential equation 



„,, 2 , sin(2/) 

f" + _ f ' !: t > = o 



(9) 



where p — r/(T — t) and ' = d/dp. For t < T we have < p < oo. 

It is sufficient to consider equation (|9|) only inside the past light cone of the point (T,0), 
that is for p € [0, 1]. The regularity of solutions at the endpoints of this interval enforces the 
following behavior 



and 



f(p) ~ ap as p -> 0, 



f(p)~- + b(l-p) as 



(10) 



(11) 



where a and b are arbitrary constants. At each endpoint the parameters a and b determine 
unique local solutions. One can show that there is a countable sequence of pairs (a n ,b n ) for 
which the corresponding solutions, denoted by f n {p), are globally regular in the sense that they 
satisfy both boundary conditions (|Io| ) and jTl| ) and are smooth for p S (0, 1). 



n 





1 


2 


3 


4 




2 


21.757413 


234.50147 


2522.0683 


27113.388 


b n 


1 


-0.305664 


0.0932163 


-0.0284312 


0.0086717 



Table 1: The parameters of solutions shown in Fig. 1. 



These solutions can be smoothly extended for p > 1 by solving (9) with the initial condition 
(11). One can show that the asymptotic behavior for p —> oo is 

/«(p)~c ri + ^+0(i), (12) 
p p 2 

where c„ — > n/2 as n — > oo. The countable family /„ was discovered numerically by Aminneborg 
and Bergstrom Q ; recently its existence was proven rigorously via a shooting argument || . The 
integer index n = 0, 1, . . . denotes the number of intersections of f n (p) with the line / = n/2 
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Figure 1: The first five self-similar solutions. 



on the interval p € [0, 1). The "ground state" solution of this family is the TS solution /q = 
2 arctan(p). The solutions f n with n > can be obtained numerically by a standard shooting-to- 
a-fitting-point technique, that is by integrating equation (9) away from the singular points p = 
and p = 1 in the opposite directions with some trial parameters a and b and then adjusting these 
parameters so that the solution joins smoothly at the fitting point. The profiles of solutions 
generated in this way (for n < 4) are shown in Fig. 1; the corresponding parameters (a„, b n ) are 
given in Table 1. 

The role of self-similar solutions /„ in the evolution depends crucially on their stability with 
respect to small perturbations. In order to analyze the linear stability of the self-similar solutions 
it is convenient to define the new time coordinate r = — ln(T — t) and rewrite Eq. (^J) in terms 
of r and p 



2 2 sin(2u) 

U TT + Ut + 2pU pT - (1 - P"){Upp + -u p ) H ^ — = 0. 

P P 



(13) 



The standard procedure is to seek solutions of ( |l3| ) in the form m(t, p) — f n (p) + w(r,p). Ne- 
glecting the 0(w 2 ) terms we obtain a linear evolution equation for the perturbation w(t, p) 



jl ^9 n a v 2 \ , 2cos(2/») 

W TT +W r +2p Wp T - (1 - p ){Wpp H Wp) H ^ w = 0. 

P P 

Substituting w(t 1 p) — e Xr v\(p)/p into (n3) we get an eigenvalue problem 



(1 - p 2 K + 2AM + A(A - l)v x + 



2cos(2/ n ) 



v x = Q. 



(14) 



(15) 



Near p = the leading behavior of solutions of ( |15[) is v\(p) ~ p a where a(a — 1) = 2, so to 
ensure regularity we require 



v\(p) 



as p — > 0. 



(16) 
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Near p = 1 the leading behavior is v\(p) ~ (1 — p) 13 where (3((3 — 1 + A) = 0. The behavior 
corresponding to [3 = 1 — A is not admissible (unless A = 1), so regular solutions must have 
(3 = 0. Then we have (up to a normalization constant) 

vx(p) ~l+ 2 + X ( 1 x ~ X) (l-p) + 0((1 - p) 2 ) as (17) 

To find the eigenvalues we need to solve Eq. (|l5|) on the interval p e [0,1] with the boundary 
conditions ( |l6| ) and ([l7]). We did this numerically (for n < 4) by shooting the solutions from 
both ends and matching the logarithmic derivative at a midpoint. Given an eigenvalue A, the 
eigenfunction v\ (p) can be extended for p > 1 by solving ( |l5| ) with the initial condition (|l7]) . Our 
numerical results strongly suggest that the solution /„ has exactly n + 1 positive eigenvalues 
(unstable modes). We denote them by (fc = 1, ... ,n + 1) where > A^ > ... > 
X„2i — 1. For example, for n ~ 1 we have A^ « 6.333625, a% = 1; for n = 2 we have 
A^ 2) Pa 59.07, A^ 2) » 6.304, A^ 2) = 1. 




Figure 2: The profiles of unstable modes around the solution /i(p). The "real" unstable mode 
(solid line) corresponds to the eigenvalue w 6.333625. The gauge mode (dotted line) has 
A2 — !• For better visualization both plots are normalized to the same slope at the origin. 



For every n the lowest positive eigenvalue A 
to the freedom of choosing the blowup time T. 
In terms of the similarity variables r = — ln(T 



= 1 corresponds to the gauge mode which is due 
. To see this, consider a solution f n (r/(T' ~ t)). 
— t) and p = r/(T — t), we have 



4t^)=^(iT^) where ( = T '- T - (18) 
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In other words, each self-similar solution f„(p) generates the orbit of solutions of ( |l3| ) f n { 1+ p eeT )■ 
It is easy to verify that the generator of this orbit 

satisfies (|l4|), thus v n = f' n (p) satisfies ( |l5| ) with A = 1. Note that this eigenfunction has exactly 
n zeros on p e (0,1) (since /„ has n extrema). For a standard Sturm-Liouville problem this 
would imply the existence of n eigenvalues above A = 1. It seems feasible to prove a similar 
result in the case of (Il5|), however we do not pursue this issue here. 



= e T pf n {p) 



(19) 



3 Static solutions 

Static solutions of Eq. (||) can be interpreted as spherically symmetric harmonic maps from the 
Euclidean space K 3 into S 3 . They satisfy the ordinary differential equation 

,, 2 , sin(2u) 

u" + -u' = °> (2°) 

r r z 

where now ' = d/dr. The obvious constant solutions of ( pfj| ) are u = and u = ir; geometrically 
these are maps into the north and the south pole of S 3 , respectively. The energy of these maps 
attains a global minimum E = 0. Another constant solution is the equator map u = it/2 but 
this solution is singular and has infinite energy. The scale invariance of ( pp| ) excludes existence 
of nontrivial regular solutions with finite energy. However, there exists a regular solution with 
infinite energy, denoted here by us(r), which behaves as 



for r — > 0, 
r 



,is(r) ~\f + ^sin(^ln. + <5) for" -- l2 " 



The existence of this solution, shown in Fig. 3, can be easily proven using x — ln(r) which 
transform (20) into a damped pendulum equation ||. 

Note that by dilation symmetry, the solution us(r) generates the orbit of static solutions 
Ug(r) = us{ar). 

We consider now the linear stability of the static solution us- Inserting u(t,r) = ug(r) + 
e lkt v(r) into (||) and linearizing, we get the eigenvalue problem 

-v" - -v' + V(r)v = k 2 v, V(r) = 2c ° 8 ^ Us \ (22) 

If the singular 1/r 2 part is subtracted from V , then (p2|) becomes the p-wave radial Schrodinger 
equation in the regular potential V reg = V(r) — 2/r 2 . This potential has infinitely many bound 
states as can be shown by the following standard argument. Consider the perturbation induced 
by the scaling transformation 

v(r) = ±u s (r) =ru' s (r). (23) 
da a=l 

This is an eigenfunction to zero eigenvalue (so called zero mode). Since us(r) has infinitely many 
extrema, the zero mode has infinitely many nodes which implies by the standard result from 
Sturm-Liouville theory that there arc infinitely many negative eigenvalues, and eo ipso infinitely 
many unstable modes around us(r). We found numerically that the "most unstable" mode has 
the eigenvalue k 2 — —0.061306. The spectrum of perturbations around the rescaled solutions 
u s (r) is obtained by scaling v(r) — > v(kr), k 2 — > a 2 k 2 . 
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Figure 3: The static solution us(r). 



To summarize, there exists the static solution us(r) (and the continuous family of its rescal- 
ings Ug(r)) which has infinite energy and infinitely many unstable modes. Can such a beast 
play any role in dynamical evolution? The answer is not clear to us. The point is that both the 
infinite energy and infinite instability of us(r) have an origin in the far-field behavior so it does 
not seem impossible that solutions Ug(r) truncated at some radius appear as local attractorsQ. 

An alternative way of looking at this issue is to consider solutions of ( p0| ) in a finite region 
r < R, that is harmonic maps from a ball B 3 (R) into S 3 . In the Dirichlet case, u(R) = c, the 
number of such solutions depends on the value of a constant c - this was discussed in detail by 
Jager and Kaul In the Neumann case, u'(R) — 0, which might be more relevant in dynamics, 
there exists a countable family of finite energy regular solutions Uk (r) . They are given by 

«*(»•)= «s*(r) with a k = ^, (24) 

where is the fc-th Neumann point of «s(r), that is a point where u' s (rk) = (fc = 1,2,...). 
By construction the solution Ufc(r) has k — 1 extrema on r G (0, R). By the same Sturm- 
Liouville theory argument as above, one can show that a truncated solution u k has exactly h — 1 
instabilities so a priori it might appear as a codimension (k — 1) local attractor in the dynamical 
evolution. 

In passing we remark that a similar structure of static solutions arising in the chemotaxis 
problem was discussed by Brenner et al. 0. 



4 Numerical results 

In this section we describe the results of our numerical simulations of the Cauchy problem (^J)-([|). 
The main goal of these simulations was to identify the generic final states of evolution (stable 

'In a recent paper |^ Liebling, Hirschmann, and Isenberg claim to have seen the solutions Ug(r) at the 
threshold for singularity formation in the evolution of very special initial data of noncompact support. We have 
serious misgivings about this result, in particular we do not understand the discussion of "critical" solutions 
which are not intermediate attractors. 
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attractors) and determine the boundaries between their basins of attraction. We emphasize 
that the convergence to attractors (which is due to radiation of energy to infinity) is always 
meant in a local sense. Before going into details, we would like to say a few words about 
the numerical techniques we employed. The simulations were performed by two different finite 
difference methods. The first method was based on an adaptive mesh refinement algorithm. This 
code allowed us to probe the structure of solutions near the singularity with good resolution. The 
second method, designed specially to study the convergence to self-similar solutions, solved the 
Cauchy problem for Eq.(|i"3|) on a fixed grid. In this case there was no need of mesh refinement 
because the convergence to self-similar profiles is a smooth process in similarity variables. The 
main difficulty of using similarity variables is that we do not know the blowup time T in advance, 
which means that we have to deal with the gauge mode instability. To suppress this instability 
(that is, to guess a blowup time T) we fine-tuned an extra parameter in the initial data. The 
fact that the two independent numerical techniques generated basically the same outputs makes 
us feel confident about our results. 

We remind that initial data of finite energy can be classified according to the topological 
degree of the wave map at a fixed time (which is a map from topological S 3 into S 3 ). Since 
the degree is preserved by evolution, the Cauchy problem breaks into infinitely many topological 
sectors. The nonzero degree data are not small by definition, and we conjecture that they always 
develop singularities. Thus, from the point of view of studying the threshold for singularity 
formation, only degree zero data are interesting so most of our discussion is focused to such 
data. A typical example is an ingoing "gaussian" 



In agreement with rigorous results of flq, |ll|] we found that if the initial data are sufficiently 
small then the solution disperses, that is it converges uniformly on any compact interval to the 
"vacuum" solution u = 0. In contrast, large initial data develop singularities in finite time - 
this manifests itself in an unbounded growth of the gradient of solution at r = 0. The precise 
character of blowup will be described below. 

Threshold behavior 

In order to determine the boundary between two generic asymptotic states of evolution, disper- 
sion and collapse, we considered the evolution of various interpolating one-parameter families of 
initial data (</>(r, p), ip(r, p)), that is such families that the corresponding solutions exist globally 
if the parameter p is small and blow up if the parameter p is large. Along each interpolating 
family there must exist a critical parameter value p* which separates these two regimes. Given 
two values p sm aii and pi ar ge, it is straightforward (in principle) to find p* by bisection. Re- 
peating this for many different interpolating families of initial data one obtains a set of critical 
data which by construction belongs to the threshold. In order to figure out the structure of the 
threshold one needs to determine the flow of critical data. The precisely critical data cannot be 
prepared numerically but in practice it is sufficient to follow the evolution of marginally critical 
data. We found that the flow of such data has a transient phase when it seems to approach the 
self-similar solution /i(r/(T — £)) for some T (see Fig. 4). This behavior is universal in the sense 
that it is independent of the family of initial data; only the parameter T depends on the data. 

This kind of behavior can be naturally explained as follows. As we showed above, the self- 
similar solution fi has exactly one unstable mode (apart from the gauge mode) - in other words 
the stable manifold of this solution, Wg(/i), has codimension one and therefore generic one 
parameter families of initial data do intersect it. The points of intersection correspond to critical 
initial data that converge asymptotically to f±. The marginally critical data, by continuity, 
initially remain close to Ws(fi) and approach f\ for intermediate times but eventually are 




(25) 
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repelled from its vicinity along the one-dimensional unstable manifold (see Fig. 5). Within this 
picture the universality of marginally critical dynamics in the intermediate asymptotics follows 
immediately from the fact that the same unstable mode dominates the evolution of all solutions. 
More precisely, the evolution of marginally critical solutions in the intermediate asymptotics can 
be approximated as 

u(t, r) = fi(p) + c(p)e Xr v(p)/p + decaying modes, (26) 

where p = r/{T — t), r = — ln(T — t), and A = A^ sw 6.3336. The small constant c(p), 
which is the only vestige of the initial data, quantifies an admixture of the unstable mode - for 
precisely critical data c(p*) = 0. The "lifetime" t* of the transient phase during which the linear 
approximation ( p6[ ) is valid is determined by the time in which the unstable mode grows to a 
finite size, that is c(p)e Xr ~ 0(1). Using c(p) w c'(p*)(p — p*), this gives r* ~ — j In \p — p*\. 




Figure 4: The plot of u(t, p) against p from the evolution (in similarity variables) of two 
marginally critical gaussian-type initial data of the form (^||), one subcritical (solid line) and 
one supercitical (dashed line). These data are identical, except for the amplitudes which differ 
by 10 -17 , so the solutions practically coincide until the last frame. The influence of the gauge 
mode instability is minimized by fine-tuning the width of the gaussian. The convergence to the 
self-similar profile f\{p) (dotted line) is clearly seen. In the last frame the two solutions depart 
from the intermediate attractor in the opposite directions. 
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Figure 5: Departure of the supercritical solution shown in Fig. 4 from the intermediate attractor. 
The T-derivative of the solution is shown to coincide (for small p) with the suitably normalized 
unstable mode around f\. 

Universality of blowup 

We address now the question: what is the shape of solutions as they approach the singularity? We 
consider first the kink-type initial data of degree one, for example u(0, r) = <f>(r) = 7rtanh(r/s). 
We found that such data always blow up in a finite time T and the asymptotic shape of solution 
near r — approaches the TS solution f (r/(T — tj) as t — > T~ . In this sense the singularity 
formation can be considered as local convergence to /o. This is shown in Figs. 6 and 7. 

We have observed the same behavior in other topological sectors, in particular in the case of 
supercritical degree zero data. In Figs. 8 and 9 we show the formation of a self-similar singularity 
in the collapse of slightly perturbed solution f\. On the basis of these numerical observations, 
we conjecture that for a large set of solutions which blow up in finite time, the asymptotic shape 
near the singularity is given by the self-similar solution / . More precisely, for such solutions 
there exists a time T such that 



inside the past light cone of the point (T, 0). 

Note that in the case of self-similar blowup the energy does not concentrate at the singularity; 
in fact the energy inside the past light cone of the point (T, 0) decreases linearly with T — t. This 
suggests that the solutions can be continued beyond the blowup time. Indeed, in Figs. 8 and 9 
we show solutions just after the blowup. At r — the solution u(t, 0) jumps from to 2w as t 
crosses T . Since u — and u — 2tt correspond geometricaly to the same point, namely the north 
pole of S 3 , the solution passing through the singularity remains smooth everywhere, except at 
one point (0,T). Moreover, the solution retains the self-similar profile at least for some time 
after the blowup. 




as t 



(27) 
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Figure 6: The evolution of kink-type initial data u(0,p) = 7rtanh(p/s) in similarity variables. 
The solution (solid line) converges to the Turok-Spergel solution / ( 1+ ^ eT ) (dotted line). By 
fine-tuning the parameter s, an admixture of the gauge mode instability quantified by e was 
made very small, e = — 0.0085 e~ 13 , so for times r < 13 the profile fo{ 1 f T ) is practically static. 




Figure 7: The same solution as in Fig. 6 at later times when the gauge mode instability shows 
up. The solution follows the moving attractor /o ( ) (dots). 
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Figure 8: The last stages (|T — t\ < 10~ 5 ) of collapse of marginally supercritical initial data 
(the solution f\ was gently "pushed" towards collapse). The arrow indicates the direction of 
increasing time. The rapidly evolving inner region and the almost frozen outer region can be 
clearly distinguished - this is a typical situation in the formation of a localized singularity. The 
numerical solution passes through the blowup in an almost continuous manner - only the point 
u(r = 0, t) jumps from to 2ir as t crosses T. 




Figure 9: Evidence of universal self-similarity of blowup. The profiles just before and after the 
blowup are shown to coincide with the Turok-Spergel solution and its reflection. 
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5 Conclusions 



We have studied the Cauchy problem for corotational wave maps from 3 + 1 Minkowski spacetime 
into the 3-sphere. We found that self-similar solutions play a special role in the dynamical 
evolution. The stable self-similar solution (Turok-Spergel solution) determines the asymptotic 
profile of solutions that blow up in finite time. The self-similar solution with one instability 
plays the role of a critical solution, that is, its stable manifold separates solutions that blow up 
from solutions that disperse. Of course, it is impossible to explore numerically the whole phase 
space, so the complete picture of singularity formation and critical behavior might be richer 
than the one sketched here. In particular our analysis leaves open the question about the role 
of a family of static solutions. Although we have not systematically investigated the nontrivial 
topological sectors of the model, we anticipate a rich phenomenology of singularity formation 
for solutions with high degree; for example we have observed such solutions evolving (in a weak 
sense) through a sequence of blowup times Ti. 

In our opinion the most interesting open question is: why the large data solutions become 
self-similar near the singularity? We think that this problem should be approached in similarity 
variables in which the problem of blowup translates into a question of asymptotic behavior 
as r — > oo. Note that the evolution equation expressed in similarity variables ( |l3| ) resembles 
the wave equation with damping. It is thus natural to seek a Lyapunov functional, that is a 
functional that decreases in time on solutions. If such a functional exists then its minima are the 
candidates for generic asymptotic states of evolution, while its saddle points are the candidates 
for positive codimension attractors. Now, the self-similar solutions /„ restricted to the interval 
p £ [0, 1] are the critical points of the functional 



Although we were unable to show that this is a Lyapunov functional for Eq. ( |l3| ) , we believe that 
the mechanism suggested here is responsible for asymptotic self-similarity of blowup. 
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